stability in MaVaN Models 
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Mass-varying neutrino (MaVaN) models propose a source of dark energy in a new scalar field 
called the acceleron. Recent work has shown that nonrelativistic neutrino fields in these theories are 
unstable to inhomogeneous fluctuations, and form structures that no longer behave as dark energy. 
One might expect that in multiple-neutrino models, the lighter species could continue to act as a 
source for the acceleron, generating dark energy without the help of heavier species. This paper 
shows that by considering the evolution of the acceleron field for a large class of models, the result 
of any component becoming unstable is that all components become unstable within a short time 
on cosmological scales. An alternate model employing a second scalar field in a hybrid potential is 
shown to have stable MaVaN dark energy even in the presence of unstable heavier components. 



INTRODUCTION 

Recent results in precision cosmology have created a 
number of puzzles. Among the unexplained phenomena 
are several apparent coincidences in which the energy 
density of two components are comparable despite having 
different redshift properties. Another is the existence of a 
negative pressure cosmological fluid that accelerates the 
rate of expansion of the universe. The model of mass- 
varying neutrinos (MaVaNs), introduced by Fardon et al. 
in [ij, suggests that neutrino and dark energy densities 
have tracked each other throughout the lifetime of the 
universe through a new scalar field called the acceleron. 
The energy density of the scalar potential of this new 
field contributes to the dark energy. The authors showed 
that this new field is capable of explaining the present 
cosmological expansion. 

Since the original MaVaN proposal, several authors 
have investigated the stability of the dark sector, com- 
posed of neutrinos and the scalar field, under perturba- 
tions to the neutrino density Both groups subject 
the model to a hydrodynamic analysis. In this approxi- 
mate picture, the speed of sound squared in the cosmo- 
logical fluid, given by Q , 
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where z is the redshift. However, for nonrelativistic neu- 



trinos, is a negative quantity while the right hand side 



dz 



is positive for w close to -1. Since at least one neutrino 
must be nonrelativistic today, the dark sector appears 
unstable. 



* Email: cspitzer@u. washington.edu 



Additionally, Afshordi, et al. perform a stability anal- 
ysis using kinetic theory to account for neutrino stream- 
ing 01 • This analysis also shows that perturbations in 
the neutrino field become unstable when the mass of the 
neutrino is of order the neutrino temperature. The ratio 
of mass to temperature when the instability occurs is a 
function of the acceleron potential, but is approximately 
7 for the potentials considered in this paper. Afshordi, et 
al., also examine the result of undergoing a phase transi- 
tion in a neutrino component, and suggest that the unsta- 
ble neutrino field may rapidly form nonlinear structures 
termed neutrino nuggets, which redshift as dark matter 
and cannot drive the cosmic expansion. 

In this paper, we will assume that the neutrinos are 
initially relativistic. As the universe expands and cools, 
the neutrinos become less relativistic until their mass and 
temperature are approximately equal. At this point, we 
assume that they decouple from the scalar field into some 
structure such as the one described by Afshordi et al. The 
resulting dark matter will not provide a large contribu- 
tion to the energy density, and we will not include the 
contribution in the dark sector energy. Note that in this 
paper dark sector will always refer to the energy density 
contributions of the stable neutrinos and the scalar field. 

In a model with multiple neutrinos, the dark energy 
may still be driven by relativistic species even after the 
heavier components have become unstable. However, the 
coupling between the neutrinos and acceleron create a 
feedback mechanism. The shift in the acceleron expecta- 
tion value when a neutrino becomes unstable can be suffi- 
cient to necessarily change the mass of another species so 
that it goes from relativistic to nonrelativistic, causing it 
to become unstable as well. This cascaded instability can 
make all neutrinos unstable at about the same time that 
the heaviest neutrino becomes unstable. This paper will 
examine a particular class of models, and will show that 
see-saw MaVaN models with flat scalar potentials suffer 
from precisely this problem. The timing between the in- 
stability in successive neutrinos is strongly dependent on 
the flatness of the scalar potential, but a flat potential is 
also required to generate the observed dark energy. This 
result may point toward models that do no suffer from 
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a cascaded instability, and this paper concludes with a 
simple example. 



SEE-SAW MODELS 

Consider a model with n active neutrinos in which each 
neutrino is paired with a sterile counterpart, which is 
coupled to a new scalar field. 
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where Vi are the active neutrinos, Ni are the sterile neu- 
trinos, and A is the acceleron scalar field. Mi and A de- 
scribe coupling strengths. If we assume that the normal 
see-saw limit holds, (XA) ^ Mi, the system reduces to 
an effective Lagrangian describing active neutrinos with 
a Majorana mass term. 
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where Mf = M^'^/A, and rrii is the effective mass. 
The energy density of the dark sector is 

Pd = Py + V{A) 



(4) 



(5) 



where V is the potential of the acceleron field. Assum- 
ing that the neutrino distribution function is a stretched 
thermal distribution, the value of the equation of state 
for this dark sector is (one derivation is given in (^) 
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where i runs over the active neutrino species, and we have 
defined the distribution function and its derivative. 
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Fardon, et al., show in that the system remains 
very close to the minimum of the effective potential and 
evolves adiabatically. For the model above, this mini- 
mization conditions becomes 
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where y — ^p-. 

Commonly employed forms for the potential include 
small power law {V = BA'',k ^ 1), logarithmic {V — 
Blog{A/Ao)) and quadratic {V = BA^). After starting 
with the small power law case, it will be easy to generalize 
to all cases by assuming dV/dA = BkA''^^ with B and 
k unrestricted. 



APPROXIMATIONS 

The expectation value of the acceleron at a particu- 
lar value of z is determined by the minimization equa- 
tion Unfortunately, this equation in general does 
not have a closed form solution. To examine the behav- 
ior of this equation it is useful to approximate the result 
in three ranges: relativistic {m^ <C T), quasirelativis- 
tic (m^ ^ T), and nonrelativistic (m^ S> T). In these 
approximation, the minimization equation becomes 
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with the unitless 0(1) integrals 
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A comparison of these approximations to unapproxi- 
mated numerical simulation for the two neutrino models 
discussed in the next section is shown in in figure^ This 
figure shows the contribution of the derivative of the neu- 
trino term to the minimization equation, which are the 
right hand sides of equations IjlOt lll() . The relativistic 
and nonrelativistic approximations (dashed curves) are a 
good match to the numerically calculated value in their 
respective regions of validity at high and low redshift. 
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FIG. 1: Tlie derivative of neutrino energy density with respect 
to A. The solid line is a numerically determined result, and 
the dashed lines are approximations in the nonrelativistic (low 
z) and relativistic (high z) regimes. 
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TWO ACTIVE NEUTRINOS, SMALL POWER 
LAW 



As a warm-up, consider the case of two active neu- 
trinos, with mi 3> m2, and a power law potential with 
< k 1. To first order, the derivative of the potential 
becomes (accurate to 10% for typical small values of k) 



dV 
'dA 



Bk 



(16) 



At large z, both neutrinos will be relativistic. As the 
universe cools, the neutrino temperature decreases while 
mass increases, and at some value of z neutrino 1 becomes 
quasirelativistic, and will become unstable. Due to the 
original mass hierarchy, neutrino 2 is still relativistic at 
this point. 

Now consider the acceleron expectation value and neu- 
trino masses both before and after the transition of neu- 
trino 1 to a dark matter phase. Before the transition, 
when neutrino 1 has mass of order the temperature, we 
have 



^1, be/ore 



T ^ A, 



before 



nT'2,before 



A 



before 



Ml 
Ml 



(17) 
(18) 



So our mass hierarchy implies Af| <C Aff . 

The minimization equation ((SJ, with the approxima- 
tions from H10() and H12() , yields the value of the acceleron 
prior to the transition, 



A 



before 



(19) 



Since the acceleron will suddenly change value when 
neutrino 1 becomes unstable and stops sourcing it, we 
do not know a priori whether neutrino 2 will be R, NR 
or QR after the transition. By trying each in turn, we 
quickly find that only the QR assumption is consistent. 
For instance, an assumption that neutrino 2 stays rela- 
tivistic yields a value for m2 that is of order the temper- 
ature, which violates the assumption. Looking carefully 
at the QR case, from the mi ~ T condition before the 
transition, we have the relationship 



(20) 



T 



which provides the temperature at the time of transition 
in terms of the parameters of the system. Solving for the 
acceleron and neutrino 2 mass after the transition, we 
have 



^after 




(21) 
(22) 



The picture that emerges from this small example is 
that after the first neutrino becomes unstable and de- 
couples from the acceleron, the acceleron assumes a value 
that pushes the second neutrino into a quasirelativistic 
region. Once this occurs, the second neutrino will also 
soon become unstable. In the end there is nothing left to 
drive the dark energy. 



GENERALIZED MODELS 

The simple two-neutrino, small-power law model gen- 
eralizes easily to include a larger number of neutrino 
species and different potentials. Additional neutrinos do 
not improve the stability picture, but decreasing the flat- 
ness of the potential does. 

First consider the addition of a third neutrino that is 
much less massive than neutrino 1. This requires Af| <C 
Mf and Af| < Aff . After the transition, the value of 
the acceleron becomes 



Aafter 



(23) 



As a result, the mass of the second and third neutrinos 
become 
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Neutrino 2 is relativistic only if M2 ^ M3, but simi- 
larly neutrino 3 is relativistic only if M3 ^ M2. Since 
we can not satisfy both these conditions, at least one of 
the remaining neutrinos is quasirelativistic, and becomes 
unstable. 

Using similar arguments, it is easy to see that for the 
general case of n neutrinos, all neutrinos will become 
unstable within a short period of each other. 

Now consider a more general potential, V = BA^ (for 
an arbitrary k), for the two- neutrino case. The mass of 
neutrino 2 after neutrino 1 becomes unstable is 



m2Mfter 



7^ Ml' 



T 



(26) 



In this case, to keep neutrino 2 relativistic, we require 

2k 

M2 \ " + 



Ml 



< 1 



(27) 



Note that this does allow the light neutrino to continue 
to drive the acceleron, but only if the acceleron potential 
is flatter than the small-fc power we examined above. Un- 
fortunately, this potential predicts a dark energy equa- 
tion of state that is in conflict with observation. The 
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value of k is connected to the present value of the equa- 
tion of state by 



Requiring a value of wq close to —0.9 yields a value of k 
that lies in the region where the lighter neutrinos become 
unstable very quickly. Conversely, values of k that are 
large enough to keep the lighter neutrinos relativistic also 
predict a value of wq too large. 

NUMERICAL SIMULATION 

We verified the above relationships using an unapprox- 
imated numerical simulation. An example of the depen- 
dence of the stability of lighter neutrinos on the steepness 
of the acceleron potential is shown in figure El The plots 
shows the evolution of the equation of state for the neu- 
trino and dark energy components, as given in JHJ, calcu- 
lated numerically, of a two neutrino system for two values 
of k. Note that this w does not include contributions from 
either the dark matter terms that resulted from unsta- 
ble neutrino components, or from components that were 
not included in the dark sector as defined above. As 
k increases, the lighter neutrino is longer-lived, and can 
continue to drive the dark energy closer to today (z — 0). 
However, this also results in w approaching a disallowed 
value. 

Comparing the set of approximations made in simpli- 
fying the minimization equation to the numerical results 
show that they are accurate to 15-50%. Noting that 
the order-of-magnitude arguments above already contain 
multiplicative 0(1) corrections, the contribution of nu- 
merical inaccuracy does not affect our conclusions. 

HYBRID MODEL 

The condition for the stability of lighter neutrino 
species, given in equation lf?fjl . suggests that we need 
to find a model that does not require a fiat potential to 
generate dark energy. There are several ways to achieve 
this, and here we will discuss a particularly simple ex- 
tension. This new model employs a second scalar field 
that is not directly coupled to the neutrinos but provides 
large contribution to the dark energy. Since the configu- 
ration of the potential is borrowed from hybrid inflation 
(see fi]), the new model is called the "Hybrid MaVaN" 
model. 

The hybrid model includes a scalar field, a, that acts 
as the waterfall field. The potential is 

V = b^A^+g^A^a^ + {h^~aa^f (29) 

where b,g,h and a are new coupling constants. The pre- 
cise details of the potential are unimportant, as long as 
it supports a false minimum as discussed below. 
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FIG. 2: Simulation of the dark sector (neutrinos and dark 
energy) equation of state of a two neutrino system under the 
power law with two values of k. The dashed line is the result if 
instability is ignored, while the solid line removes an unstable 
component from the dark sector. The solid lines end when 
all neutrinos have become unstable. The contributions of the 
dark matter formed by unstable neutrino components is not 
included here. 

The hybrid potential has two distinct regions of be- 
havior under variation of a. If 2ah'^ > g^A'^, then there 

are two minima at ± -^/^i^E^ilZE^ Otherwise, there is a 
single "false minimum" at 0. Forcing the field into this 
false minimum by requiring 

A > V2^h/g (30) 

the potential becomes 

V -> ^2^2 + (31) 

which includes a cosmological-constant type term h'^. 
This term dominates the acceleron contribution to the 
potential if h:^ V^b/g, and from equations (0) and © 
may also dominate over the neutrino contribution to the 
energy density. If we also assume the see-saw condition, 
A 3> rrii/X, then the model described in previous sec- 
tions can be used without any change other than using 
the form of the potential in equation l|31|) . 

The quadratic dependence on A in equation H31|l means 
that the stability condition for the lighter neutrinos in 
equation (|27|l is easily satisfied. The allowed parameter 
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range is quite large, and it is straightforward to find co- 
efficients that are produce observationally allowed values 
of neutrino mass and equation of state. The numerical 
simulation of the evolution of one such model, with a hi- 
erarchy of masses and h = 0.06 is shown in figure^ The 
lighter two neutrinos stay relativistic until z = 0, despite 
the instability in the massive neutrino. Note this model 
achieves both the stability of the lighter neutrinos and 
has a dark sector equation of state of w = — 1 at z near 
0. 

^ Hybrid 




FIG. 3: Simulation of the equation of state of a three neutrino 
hybrid potential system. The kink at z = 200 occurs when 
the heaviest neutrino becomes unstable. The remaining two 
neutrinos continue to be relativistic and stable until z—Q. 



There are a number of possible resolutions. One is to 
increase the curvature of the scalar potential. In models 
with a single scalar field, this makes it difficult for the 
scalar field potential to form dark energy. However, this 
is easily remedied by including a second scalar field. A 
simple example is illustrated above in the hybrid MaVaN 
model. A similar potential arises naturaljv in supersym- 
metric models, such as those in and 

Another solution is to modify the theory so that the 
dark sector never reaches a state where the adiabatic con- 
dition applies. Such models do not suffer from the in- 
stability described above. One such theory is presented 
in 9J. 

Reducing the dependence of the acceleron on the heavy 
neutrino components also forms a class of possible so- 
lutions. If the acceleron is decoupled from each heavy 
component before it becomes unstable, the acceleron ex- 
pectation value is not quickly driven to a new scale. The 
mass of the lighter neutrinos would be largely unaffected, 
avoiding instability. 
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CONCLUSION AND DISCUSSION 

Several authors have found that a nonrelativistic Ma- 
VaN neutrino field is unstable to inhomogeneous fluc- 
tuations. By considering the evolution of the acceleron 
expectation value when a neutrino field becomes unsta- 
ble, we have argued above that all neutrino fields in the 
theory are susceptible to a cascaded instability in which 
they all become unstable at nearly the same time. This 
occurs as long as the scalar potential has a nearly flat 
dependence on the acceleron, and the neutrino masses 
vary inversely with the acceleron. Including a very light 
neutrino is not sufficient to avoid this problem. Since 
there are at least three neutrino species, and the atmo- 
spheric neutrino deficit requires at least one mass scale 
above leV, the instability poses a constraint on all phys- 
ical MaVaN models. 
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